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ABSTRACT 



Aims. We study galaxy clustering and explore the dependence of galaxy properties on the the environment up to a redshift z ~1, on the basis of 
a deep multi-band survey in the Chandra Deep Field South. 

Methods. We have developed a new method which combines galaxy angular positions and photometric redshifts to estimate the local galaxy 
\£) number-density. This allows both the detection of overdensities in the galaxy distribution and the study of the properties of the galaxy popula- 
f"^ ■ tion as a function of the environmental density. 

Results. We detect two moderate overdensities at z~0.7 and z~ 1 previously identified spectroscopically. We find that the fraction of red galax- 
ies within each structure increases with volume density, extending to z~l previous results. We measure "red sequence" slopes consistent with 
the values found in X-ray selected clusters, supporting the notion that the mass-metallicity relation hold constant up to z~l. 
O ' Conclusions. Our method based on photometric redshifts allows to extend structure detection and density estimates up to the limits of photo- 
metric surveys, i.e. considerably deeper than spectroscopic surveys. Since X-ray cluster detection at high redshift is presently limited to massive 



Or 



relaxed structures, galaxy volume density based on photometric redshift appears as a valuable tool in the study of galaxy evolution. 



Key words. Galaxies: clusters: general - Galaxies: evolution - Galaxies: formation - Galaxies: distances and redshifts - (Cosmology:) large- 
scale structure of Universe 

1 . INTRODUCTION the SDSS survey (York et al. EOOOl so that complete spectro- 

scopic galaxy redshift surveys of large areas of the sky are not 



In recent years deep multi-band galaxy surveys, together with 
the development of reliable population syntheses, and mod- 
els for galaxy formation/evolution which incorporate feedback 
from the galactic environment, have determined a dramatic in- 
crease of our understanding of the observed average evolution 



available at the boundary of the visible Universe. To minimize 

the effect of the background/foreground objects in the detection 

of galaxy overdensities, two-dimensional (2D) surveys, which 

are deeper, require additional a priori assumptions on either 

galaxy luminosity function (LF), as in the Matched Filter algo- 
of the galaxy population in cosmic time. On the other end, cms- ... /ri \ , | lr<r ,, i „ , 

, _ , . , • , • nthm (Postman et al. 1996) or the presence of a red sequence 

ters and groups of galaxies are the ideal laboratories where it ,, . , | . nni ii r.- j ji i 

., f , , . , m ... . (Gladders et al. 1998). Biases produced by these assumptions 

is possible to study the environmental effects, which cause the 



differences of the distribution of galactic types from the aver- 
age (Carlberg et al. 120011 Gomez et al. 120031 1. For this reason 



can hardly be evaluated for at high redshift. Multi-band pho- 
tometric surveys are able to provide redshifts for large statis- 
tical samples of galaxies at the deepest observational limits, 

finding overdense regions at high redshift plays a central role in , ... . . „ , ., 

° o o r j though with a lower accuracy (cr- a 0.05) as compared with 

understanding galaxy formation and evolution. Unbiased sam- " , „ „„ lw „ . . hn »»i „ .. . . 

bb J spectroscopic ones (o\ x;0.001) (Fontana et al. 2000, Bell et al. 

pies of clusters and groups of galaxies can be obtained from 3D . . ... . . . , . . . , 

* o ±- o 2004). Various algorithms exist to evaluate photometric red- 

spectroscopic surveys, where the ability of a Friend of Friend . ,. , 

fL^^ , , „, , „ ^ „ rrTTTTTn r , , • • shirts. They can rely on the position in a multi-dimensional 

(FOF) method (Huchra & Geller 1982) to find over-densities . . . „, , , .. 

1 ' colour space respect to galaxies of known redshift (Connolly 

is only limited by the uncertainty caused by the galaxy veloc- . , | lrir , g i . . . . .. . ... . ~. , , ct ,t, 

. J , i ■ ■ et. al. 1995). Alternatively it is possible to fit the observed SED 

lty dispersion. However, spectroscopic surveys reach limits of . . . /T , . . . 

, , ...... , . . with either empirical (Lanzetta et al. 1996) or theoretical tem- 

about 5 magnitude brighter than photometric ones, as e.g. in . „ ^ . i .^^o. „, ^ . , 

° plates (Giallongo et al. 1998). Photometric redshifts have been 
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2005 ). However, the present accuracy of photometric redshifts 
leaves no room for standard FOF methods. In fact, the linking 
distance dp of must be larger than the distance uncertainty 6 p i lot 
associated with the photometric redshift, otherwise even close 
neighbours, which on average appear at a r.m.s. distance of the 
order of 6 p h ot , fail to be linked. On the other hand, 6fof must 
be smaller than the average galaxy-galaxy distance d g , to avoid 
a percolation of the "friendship" path through the galaxy distri- 
bution leading to huge unphysical chains of galaxies linked by 
the FOF algorithm, even in the absence of any group or clus- 
ter structure: thus the condition 6 p i, ot < Sfof < d g must hold. 
However, the minimum uncertainty attained by photometric 
redshifts is Az ~ 0.05 • (1 + z), corresponding to S p hot >180 Mpc 
for z > 1, while the average number density of bright galaxies 
is ~ 10~ 2 M pc~ 3 , corresponding to an average galaxy-galaxy 
distance as small as d ~ 5 Mpc (Bahcall 1988), i.e. we have 
Sphot » d g . This has been discussed in detail by Botzler et 
al. [2004 who proposed an extended friends-of-friends method 
(EXT-FOF) which applies a FOF method to angular distances, 
in redshift slices which are defined on the basis of photometric 
redshifts, taking into account their intrinsic uncertainties. In the 
present paper we present an alternative approach to the use of 
photometric redshifts. Our algorithm evaluates the 3D galaxy 
density using angular positions and photometric redshifts with 
the purpose of: i) detecting galaxy over-densities in three di- 
mensions, and ii) assigning to each galaxy a measure of the 
environmental density, to extend the analysis of the environ- 
mental effects on galaxy evolution to the limits of photometric 
surveys. The method is applied to a deep photometric survey 
of the Chandra Deep Field South (Cimatti et al. 2002a ), where 
we re-discover some previously known overdensities and we 
find a clear relation between the fraction of blue and red galac- 
tic types and the local density. The paper is organised as fol- 
lows. In section 2 we present a new method to define a galaxy- 
volume density using angular positions and photometric red- 
shifts. In section 3 we apply the algorithm to a deep photo- 
metric survey of a portion of the Chandra Deep Field South 
to identify some high-z clusters and analyse the dependence of 
their galaxy population on local density. In section 4 we study 
the evolution in cosmic time of the red sequence. In section 5 
we check the reliability of the results and possible extrapola- 
tion to deeper surveys by applying the algorithm to a simulated 
galaxy catalogue. In section 6 we discuss the results and section 
7 provides a summary. Throughout the paper we adopt consen- 
sus cosmological parameters 0.^=0.1, Qm=0.3, Hq=70 km s _1 
Mpc" 1 . 

2. (2+1)-DIMENSIONAL DENSITY FROM 
PHOTOMETRIC REDSHIFTS 

As far as background/foreground effects are acceptable, namely 
at relatively low redshifts, the surface density is sufficient to 
detect galaxy clusters or to study their morphology, as in the 
case of Abell ( 1958 ) or Zwicky (Zwicky et al. |1961| ) clusters or 
in the later works of the Edinburgh-Durham Cluster Catalogue 
(Lumsden et al. 1992). A surface density 

£„ = nl{nD 2 pn ) (1) 



of galaxies can be computed considering the projected dis- 
tances D p „ to the n-th nearest neighbour, and can be used to 
study the effect of the environment on the distribution of galac- 
tic types, as done by Dressier (Dressier 1980 Dressier et al. 
I1997I ). Similarly, when spectroscopic redshifts are available, is 
possible to compute the three-dimensional (non projected) dis- 
tance D„ to the n-th nearest neighbour and derive a volume 
density 

p„ = 3n/(4nD 3 n ). (2) 

The method we propose for deep photometric surveys, where 
spectroscopic galaxy redshifts are not available, consists in 
combining, in the most effective way, the angular position with 
the (much less accurate) distance as computed from the photo- 
metric redshift. In principle it is possible to compute for each 
galaxy the distance from its neighbours from the angular po- 
sitions and photometric redshifts, once a cosmological model 
is assumed. In practise we prefer to divide the survey volume 
in cells whose extension in different directions (Aa, A6, Az) de- 
pends on the relevant positional accuracy and thus are elon- 
gated in the radial direction. Then for each point in space (i.e. 
for each cell) we count neighbouring objects of increasing dis- 
tance, until a number n of objects is reached. We define the 
density associated to the cell as 

P = n/V„ 

where V n is the volume which includes the n nearest neigh- 
bours. The choice of n is a trade off between spatial resolution 
and signal-to-noise ratio. In the following the number n is cho- 
sen in such a way that n objects are present in a single cell in 
the regions of maximum density. Notice that, in principle, the 
way of searching neighbouring objects of increasing distance 
in not univocal, since distance steps (cell sizes) in different di- 
rections can be chosen arbitrarily. The result will be a different 
smoothing, and resolution, in different directions. The choice 
we adopt is to keep the maximum resolution in transversal and 
radial directions allowed by the data. Since the uncertainty in 
the radial direction is much larger, this will correspond to elon- 
gated cells and lower radial resolution. Steps Az in the radial di- 
rection smaller than the photometric redshift uncertainty would 
uselessly increase the computing time. In counting galaxies we 
must take into account the increase of limiting luminosity with 
increasing redshift for a given flux limit. If m/,,„ is the limiting 
(apparent) magnitude in a fixed observing band, at each red- 
shift z we detect only objects brighter than an absolute magni- 
tude M/, m (z), decreasing (brightening) with z . We can assume a 
reference redshift z c below which we detect all objects brighter 
than the relevant M c = M/„„(z c ), below which we neglect the 
incompleteness. At z > z c the fraction of detected objects is: 

Mum 0(M)dM 

£* c ®(M)dM 

where <t>(M) is the galaxy luminosity function. Thus, in eval- 
uating the galaxies number density, we apply a limiting mag- 
nitude correction by assigning a weight w(z) = l/s(z) to each 
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detected galaxy of redshift z. In this way we correct the sys- 
tematic underestimate of density, caused to the increasing frac- 
tion of galaxies which fall below the brightness limit for in- 
creasing redshift. Of course the noise in the derived density 
increases as the square root of w(z), but the advantage is to ob- 
tain a density scale independent of redshift, at least to a first 
approximation. We stress that, in general, the distance modu- 
lus m - M depends not only on the luminosity distance Dl(z), 
and thus on the adopted cosmological model, but also on the 
k- and evolutionary- corrections which, in turn, depend on the 
galactic type. Moreover the luminosity function itself depends 
on both the wavelength A and cosmic time. However, as far as 
s(z) is not much smaller than one, it is possible to adopt simple 
representations of these effects, and any correction to the lim- 
iting magnitude correction will represent a second order effect. 
Practically the choice of the cell sizes is determined by the ac- 
curacy of the photometric redshifts. In fact, once the cells have 
transversal sizes which are much smaller than the radial one 
(say 1/1000), a further increase of the transversal resolution 
does not justify the corresponding increase of the computing 
time. At least, this is true as far as we do not expect structures 
which are physically strongly elongated in the radial direction. 
In the application presented in the following sections, the rel- 
ative redshift uncertainty is A = \z spe c - Zpfa»l/(1 + z) ~ 0.05 
(Cimatti et al. 2002b), corresponding to a length scale l z « 200, 
Mpc at z ~ 1. Cells sizes are arbitrary and are chosen small 
enough not to degrade the spatial resolution, while avoiding a 
useless increase of the computing time. We have adopted ~ 60 
Mpc in the radial direction and ~ 60 kpc at z ~ 1 in transversal 
direction, corresponding to Aa = A6 — 3.68 arcsec. Galaxy 
clusters or, more generally, galaxy overdensities are defined 
as connected 3-dimensional regions with density exceeding a 
fixed threshold. Once overdensities are identified, it is possible 
to analyse the properties of galaxies as a function of the local 
density. 

3. GALAXY CLUSTERS IN THE CHANDRA DEEP 
FIELD SOUTH 

The dataset used to identify distant clusters and to study the 
properties of their member galaxies is the deep photometric 
catalogue of the K20 survey (Cimatti et al. I2002al [2002bb con- 
taining photometry in the UBVRIZJK bands of a 6.38x6.13 
arcmin field in Chandra Deep Field South (CDFS) (Giacconi 
et al. 2002). The sample is limited to Iab < 25, while the pho- 
tometric depth in the other bands allows to assume that virtu- 
ally all galaxies in the catalogue have 8-band photometry, ex- 
cept a few objects with very extreme colours. We have added 
to the spectroscopic redshifts of the K20 all public spectro- 
scopic redshifts in our field from GOODS-MUSIC catalogue 
(Grazian et al. 2006 and refs. therein). The catalogue contains 
1749 galaxies among which 292 have spectroscopic redshifts 
and the remaining have only photometric redshifts. The proce- 
dure for deriving photometric redshifts and test their accuracy 
are described in Cimatti et al. I2002bl where it is shown that the 
distribution of the fractional error A = (z S pec ~ z p i w t)/(l + z spe c) 
is not gaussian. After excluding "outlayers" with A > 0.15, 
which represent less than 9% of the total, we obtain <x A = 0.05. 




spec 



Fig. 1. Photometric redshifts z p i w , versus spectroscopic red- 
shifts z spec f° r a H the galaxies in the catalogue with spectro- 
scopic observations (Grazian et al. [2006 and refs. therein). 
Dashed lines indicate the r.m.s. uncertainty 0.05(l+z). 

Table 1. Detected structures. 



# 


<*2000 


^2000 


z 


1 


03 32 20.08 


-27 47 07.20 


0.70 


2 


03 32 16.74 


-27 45 52.99 


1.00 


3 


03 32 16.51 


-27 46 45.95 


1.55 



Figure 1 shows photometric redshifts versus spectroscopic 
redshifts, with the uncertainty 0.05(1 + z) indicated by the 
dashed lines. 

We have constructed (2+l)D maps of the volume density 
p„, with n = 10, computing s(z) in equation (3) on the basis of 
a cosmologically evolving luminosity function. This has been 
taken from Poli et al. 120031 who derived the rest-frame B lumi- 
nosity function which is directly sampled until the rest-frame 
blue is observed in the K band, namely up to a redshift of about 
3.5. In their analysis Poli et al. 2003 find little density evolu- 
tion at the faint end with respect to the local values, while at 
the bright end a brightening increasing with redshift is appar- 
ent with respect to the local LF. 

The choice n — 10 corresponds to the maximum number of ob- 
jects in a single cell at high density as discussed in Section 2. 

To see how the resolution of photometric redshifts com- 
pares with real objects distribution, we show in Figure|2]the his- 
togram of photometric redshifts in the field. Two main clumps 
appear about redshifts 0.70, 1.00. A comparison with the dis- 
tribution of spectroscopic redshifts of the CDFS (see Gilli et al. 
2003 fig. 1) clearly indicates the reality of the two clumps, al- 
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Fig. 2. a) The distribution of photometric redshifts of the sam- 
ple, b) Average p\o density on the entire field, in redshift bins, 
versus z as determined by the (2D+1) algorithm using for all 
sources: i) photometric redshifts (dashed line); ii) photometric 
redshifts, or spectroscopic redshifts whenever available (con- 
tinuous line). 




3" 32 m 30 s 25 20 15 10 30 25 20 15 10 
R.A. R.A. 

Fig. 3. Isodensity contours of the surface density Sio as com- 
puted in the redshift slice 0.70 < z < 0.75, which includes the 
first detected structure (left panel), b) Same plot for the second 
structure in the redshift slice 0.90 < z < 1 . 10 (right panel). 

though the two peaks at z=0.67 and z=0.73 found by Gilli et al. 
2003 in the distribution of spectroscopic redshifts are barely re- 
solved in our photometric redshift distribution. On the basis of 
spectroscopic redshifts, Gilli et al. 12003 1 found that both struc- 
tures are spread over the field, although the latter includes a cD 
galaxy, suggesting a dynamically relaxed status (see Sarazin 
[1988l l. 

To check to what extent our results depend on the presence 
of spectroscopic redshift in our catalogue, we show in Figure 
I2J5 the average density in the field as a function of redshift, as 
computed from photometric redshifts only, or including also 
spectroscopic redshifts whenever the latter are available. The 
two curves look similar. In fact using z sp ec instead of z p hot does 
not change significantly the density, once the cell size has been 
chosen on the basis of the (much lower) z p hot accuracy. This 
means that our results do not rely on the availability of sev- 
eral spectroscopic redshifts in the field. Notice that the density 
reported in Figure [2J5, which is averaged over the entire field, 
is not used to detect structures. For this purpose we select in- 
dividual cells with density above a given threshold, and then 
we look for connected volumes. The choice of the threshold 
is an arbitrary trade off between completeness and reliability. 
From our numerical simulations (see section 5) we found that 



a threshold corresponding to about three to five times the av- 
erage density can identify richness zero Abell clusters. More 
complex simulations would be necessary to reliably evaluate 
the degree of contamination as a function of the richness class 
and redshift. In real data, from the distribution of the densities 
of individual cells we found that thresholds of 0.078 Mpc" 3 
or 0.13 Mpc~ 3 (i.e. 3 or 5 times the average) isolate 2.2% or 
0.5% respectively of the total cell number. To avoid excessive 
contamination from random density fluctuations we adopted a 
threshold pf™ sh =0.08 Mpc~ 3 which selects about 2.0% of the 
cells. In this way we identified three over-densities listed in ta- 
ble 1. We adopted a redshift independent threshold to detect 
structures of comparable density at any redshift. In principle, 
this implies an higher probability of contamination at higher 
redshifts: we plan to explore this issue in future work extend- 
ing to higher redhifts. 

The first structure, at z=0.70 is approximately centred on 
the cD galaxy, whose position, in turn, corresponds to the cen- 
tre of the extended X-ray source CDFS566 (Giacconi et al. 
2002). As noted above, we cannot resolve, in redshift space, 
the wall-like structure at z=0.67 described by Gilli et al. 2003 
which then contaminates the structure at z=0.70. In spite of 
that we find a relatively concentrated structure with full width 
at half maximum of about 0.12 Mpc. After the statistical sub- 
traction of the background/foreground field galaxies the num- 
ber of galaxies within the Abell radius Ra is N c =182, of which 
38 are between mi and mj, +2, corresponding to a richness class 
0. From the spectroscopic redshifts we can evaluate a velocity 
dispersion along the line of sight, cr p = 334 + 31 Km s _1 , for 
the 39 galaxies belonging to the peak in the redshift distribu- 
tion centred at at z=0.73, the uncertainty being computed by 

a bootstrap method. The relevant virial mass M = I is 
M = 1.19- 1O 14 M , where R PV = N(N- 1)/ Z^jRj/ is the pro- 
jected virial radius and Rjj are the projected distances between 
each pair of the N=39 galaxies (Heisler et al. 1985 Girardi 
et al. |19981 l. For the sole purpose of displaying the morphol- 
ogy of the density field we show, in Figure [3^, the isolines of 
the surface density £10 (see section 2), evaluated in the red- 
shift slice 0.70< Zphot < 0.75. A similar plot for the overden- 
sity at z^l.00 is shown in Figure |3j) where the galaxy num- 
ber within the Abell radius, with < m < 1113 + 2, barely 
reaches the formal threshold of 30 corresponding to richness 
class 0, depending on the exact location of the adopted centre. 
The spectroscopic redshift distribution suggests the existence 
of two distinct peaks at 0.97 and 1.04: the former associated 
with galaxies around the main overdensity and the latter corre- 
sponding to the south east extension. The analysis of possible 
substructures of this overdensity requires further spectroscopic 
data. The third clump we find at z — 1 .55 does not appear in 
the distribution of spectroscopic galaxy redshifts which is lim- 
ited to brighter fluxes respect to our photometric data. On the 
other hand,the accuracy of our photometric redshifts is statis- 
tically checked against spectroscopic ones only for z spec ~ 1, 
so that further data would be needed to assess the reality of 
this structure. However a peak in the distribution of the X-ray 
selected AGNs in the field is present at about z=1.55. Thus, 
as far as we can assume that distribution of AGNs traces the 
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Fig. 4. Galaxy colour distributions: at high density (pio > 
0.08Mpc~ 3 ) (shaded histogram), at low density (pio < 
0.03 Mpc~\ for the two structures at z~0.70 (left) and z~1.0 
(right). 




0.15 0.1 0.05 
p(gal/Mpc 3 ) 



0.15 0.1 0.05 
p(gal/Mpc 3 ) 



Fig. 5. Fraction of galaxies with rest frame colour B-R> 1 .25 
as a function of the volume density pio, for the two structures at 
z~0.70(left) and z~1.0 (right). Error bars represent Poissonian 
fluctuation. 



the constant rest frame colour B-R=1.25 as a boundary be- 
tween the two populations. The excess of red galaxies in over- 
dense regions respect to the field, clearly appears in Figure 
|U According to a Kolmogorov-Smirnov test, the probability 
of the null hypothesis that the colours inside and outside the 
overdensities are randomly drawn from the same distribution is 
4.2xl0~ 7 for the cluster at z=0.7 and 0.058 for the cluster at 
z=1.0. Due to the insufficient statistics, a similar colour segre- 
gation cannot be detected in the case of the z~ 1 .55 overdensity, 
which will be studied as soon as deeper photometric data will 
become available (Trevese et al. 2006). In the two former over- 
densities we can study the fraction of red galaxies as a function 
of the density pio. The result is shown in Figure|5] where both 
clusters show a decrease of red fraction as a function of pio, 
marginally significant at z=1.0 and more evident at z=0.7. This 
colour segregation is clearly related to the morphological seg- 
regation first found by Dressier 1980 for local clusters succes- 
sively extended to z=0.5 by Dressier et al. 119971 Our analysis of 
colour segregation allows to extend a quantitative investigation 
of environmental effects up to redshifts where morphological 
studies become unfeasible. 

4. COSMOLOGICAL EVOLUTION OF THE RED 
SEQUENCE 

We first study the cosmological evolution of the average red 
sequence of the field. Following Bell et al. 2004 it is possible 
to define for each galaxy a colour index C reduced to My m = 
-20 by shifting each galaxy in the C-M diagram to My resl =-20 
along the red sequence: 



C = C + a RS ■ (M 



Vrest 



20) 



large scale distribution of matter we can say that we are detect- 
ing a structure not previously seen in galaxy spectroscopic sur- 
veys. The Abell richness of the third structure at z~1.55 cannot 
be evaluated, since + 2 falls below the limiting magnitude 
m/=25. The peak in the photometric redshift distribution con- 
tains 57 objects spread along a moderate over-density cross- 
ing the field from north-west to south-east likely related to the 
above discussed large-scale structure traced by X-ray selected 
AGN (Gilli et al. 2003). The association of an environmental 
density with individual galaxies allows both a further assess- 
ment of the nature of the detected overdensities and the analy- 
sis of the relation between galaxy spectral energy distribution 
and the environment. 

A strong colour bi-modality of the colour distribution has 
recently been confirmed on the basis of a large galaxy sample 
of about 150,000 objects from the Sloan Digital Sky Survey 
(Strateva et al. 120011 and refs. therein). This bi-modality has 
been shown to maintain up to z^2-3 (Giallongo et al. 2005), 
with a local minimum in the colour distribution which evolves 
in redshift and represents the natural separation between the 
"blue" and "red" galaxy populations, the latter defining an av- 
erage red sequence of the field. Instead, Figure H shows rest- 
frame B-R colour distribution in overdense (pio > 0.08 Mpc~ 3 ) 
and underdense (pio <0.03 Mpc~ 3 ) regions both at z~0.7 and 
z~l. Following Carlberg et al. 120011 we conservatively adopt 



where C is the rest-frame colour, (U-V) resf in our case, and 
curs = d(U - V) resr /dMv rea is the slope of the red sequence, 
In practise, for a direct comparison with Bell et al. 2004 who 
adopt a different cosmology, we assume as reference magni- 
tude M Vmt = -20.7 instead of -20. A distribution of the C', 
instead of C, colours allows a better identification of the "Early 
type", or red galaxies, population which defines the average 
red sequence. Following Bell et al. 120041 we assume a constant 
slope d(U - V) res ,/dMv lra = -0.08 which is derived from a 
sample of nearby clusters (Bower et al. H992t . This assumption 
is justified by the analysis of Blakeslee et al. [2003 who find a 
constant slope for different galaxy clusters up to z=1.2 (see the 
next paragraph). At higher redshifts Giallongo et al. 2005]find 
a decrease of the red sequence slope d(U — V) rest ldMB ml from 
-0.098 to -0.062, in the two wide redshift bins 0.4-1 and 1 .3-3.5 
respectively. In the present analysis, which is limited to z<1.7, 
we neglect this change of slope and we estimate the red se- 
quence colours in three relatively narrow redshift bins. Figure 
|6](a-c) shows the rest-frame colour-magnitude diagrams for all 
galaxies in the field, selected in the three slices of photometric 
redshift, 0.7 < z < 0.8, 0.9 < z < 1.1, 1.4 < z < 1.7 and 
the relevant C' distributions. From the C' colour distribution 
we select the members of the red population as those galax- 
ies lying redwards of the natural minimum which separates the 
blue and red populations. A fit in the C-M diagram of the red 



6 D. Trevese et al.: large scale structure in the CDFS 



0.15 0.1 0.05 



1 - 




Fig. 6. Left panels: histograms of the C colour, defined in 
the text, for all the objects in the intervals 0.7 < z < 0.8, 
0.9 < z < 1.1, 1.4 < z < 1.7 (from the top). Dashed lines, 
corresponding to a local minimum, define the red and the blue 
populations. Right panels: rest-frame (U-V) vs. My diagrams: 
the continuous line represents the fit to the points belong- 
ing to the red population, with a fixed slope or ss =0.08. The 
dashed vertical line represents the reference absolute magni- 
tude My = -20.7. The dotted line indicates the valley separat- 
ing the two populations. 



population with a straight line of fixed slope a^s defines the 
colour C' RS of the average red sequence in different bins of red- 
shift. Figure 13 adapted from Bell et al. 120041 shows the colour- 
ed of the average red sequence in our field in the above red- 
shift intervals , compared with the colour of the average red 
sequence as a function of redshift obtained from COMBO 17 



data, after a correction on AC = C 



combo 



- C 

^CDFS 



= -0.08 



as derived from a comparison of the colours of those high-z 
galaxies which are in common in the 2 catalogues COMBO 17 
and K20. Our error bars represent the r.m.s. dispersion of C 
in each redshift bin. We can conclude that the AC vs. z re- 
lation in our analysis is consistent with Bell et al. 2004 up to 
z^l, moreover the point at z=1.5 lies on the linear extrapola- 
tion of their points. As a further step, we measured the slope 
of the red sequence in the two clumps at redshifts 0.7 and 1 .0. 
This has not be done for the third clump at z=1.55, due to in- 
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Fig. 7. The colour C' RS at My mt --2Q of the average red se- 
quence, as computed from COMBO 17 survey (circles), and re- 
sulting from the present analysis (triangles). The straight line 
represents a linear fit on all the data analysed in the COMBO 17 
survey (adapted from Bell et al. 120041 1 . 




Fig. 8. The slope of the rest-frame (U-B) vs. Mb in the galaxy 
clusters collected by Blakeslee et al. 2003 (open squares) and 
in the two structures detected in the CDFS (filled squares). The 
dotted line represents the average value derived by Blakeslee et 
al. 120031 

sufficient statistics. We selected galaxies with an environmen- 
tal density above a threshold pio=0.08 gal/Mpc 3 (see section 
3). We then isolated the red population from the histograms of 
the C colour which are clearly bimodal. Finally, we evaluated 
the slope of the colour-magnitude relation for this population 
in the U-B vs. B diagram, as done by Blakeslee et al. 2003 
Our results show that the slope of the red sequence is consis- 
tent with being the same in local clusters and in the two main 
overdensities (z~0.7 and z~1.0). Blakeslee et al. 2003 found 
little or no evidence of evolution of this slope, out to z=1.2. 
Our results (Fig. [8]) lie within Icr from their average value 
| < d(U - B) rest /dM Bml > | = 0.032. According to the standard 
interpretation (Arimoto & Yoshii 1987, Kauffmann & Chariot 
1998) this implies that the mass-metallicity relation holds the 
same from z=0 up to at least z=l 

5. A CHECK OF THE ALGORITHM ON SIMULATED 
CLUSTERS 

To check the reliability of the algorithm in detecting clusters of 
various types and redshifts we created a series of mock cata- 
logues including field galaxies and clusters. Field objects were 
uniformly distributed on a square of 6 x 6 arcmin centred on the 
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Table 2. LF parameters. 



m ]im =25 m lim =27 



type 


<P"o 


7 


M* 


6 


a 


Ellipticals 
Spirals 


0.0106 
0.0042 


-2.23 
-0.52 


-20.03 
-20.11 


2.72 
2.35 


-0.46 
-1.38 



cluster, with a space density and absolute magnitude distribu- 
tions assigned according to the redshift-dependent Schechter- 
like (Schechter [T976i l LF derived by Giallongo et al. 120051 



<S>(M,z)dM = O.4MO0* (z)[10' 



,0.4(M*fe)-M)]l+« ejC p|-. 



10' 



0.4(M*(z)-M) 



with 4>* = 0g(l + z) y and M* = M ( * - 5log{\ + z). Clusters are 
represented by a number-density distribution 

«(r) = « [l+Wr c ) 2 r 3/2 



(see Sarazin 1988) with a typical core radius r c = Q.25Mpc. To 
take into account the uncertainty on photometric redshifts, we 
assigned to each galaxy a random redshift with a mean corre- 
sponding to the cluster redshift and a dispersion cr : = 0.05 cor- 
responding to a velocity dispersion along the line of sight cr v = 
15000/(1 + z) Km s 1 (Hogg [2000l Once position and redshift 
are assigned to a galaxy, an absolute blue magnitude Mb is ex- 
tracted from a Schechter LF, with parameters M* B = -20.04 and 
a = -1.05 for ellipticals and M* = -19.48 and a = -1.23 for 
spirals, derived by De Propris et al. [2003 from the analysis of 
2dF clusters. The fraction of different galaxy types was chosen 
according to Goto et al. |2004[ Different values of the central 
density n are chosen to produce clusters of various richness 
classes. Table 2 reports the values of the parameters for ellipti- 
cal and spiral galaxies in the field. To evaluate the total number 
of cluster and field objects N = ®(M,z)dM we com- 

puted the limiting rest-frame magnitude in the B band given 
by: 



Mr 



= mi + C-K-E- 25.0 - 5.0 • log m D L 



mi = 25 is the limiting AB magnitude at the effective wave- 
length of the I band, K and E are, respectively, the type de- 
pendent K-correction and evolutionary correction (Poggianti 
1 19971 1, C is the relevant rest-frame (B-I) colour and Di is the lu- 
minosity distance. We assumed a fixed proportion of Sa (70%) 
and Sc (30%)in the field catalogue. Various mock catalogues 
were generated at each cluster redshifts 0.5, 1.0, 1.5, 2.0 cor- 
responding to Abell (1958) richness classes 0,1,2,3. Figure 
[9] shows the galaxies belonging to a richness class cluster, 
above a pin density threshold five times the average, as seen at 
different redshifts and limiting magnitudes. Real cluster mem- 
bers and interlopers are represented by filled dots and crosses 
respectively. For a limiting magnitude m/, m =25, the cluster is 
detected with an acceptable contamination up to z=l, while 
z=1.5 may be assumed as a detection limit. However, in the 
case of a deeper survey, with m; lm =27 comparable with the new 
GOODS-MUSIC catalogue (Grazian et al. |2006t , whose anal- 
ysis is in progress, the same richness cluster is well detected 
up to z=1.5 and still visible at z=2. 
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Fig. 9. Galaxies above an environmental density pin threshold 
5 times the average, for redshifts 0.5, 1.0, 1.5 and 2.0 and for 
catalogues with two different limiting magnitudes. Filled dots 
represent real members of the simulated richness class clus- 
ters, while crosses represent interlopers. The contamination is 
reported in each panel. 



Clearly the ability of the algorithm in separating the over- 
densities increases with their richness and their angular dis- 
tance. Thus, to perform a conservative evaluation of the algo- 
rithm, we simulated two relatively low density structures per- 
fectly aligned along the line of sight. We have considered pairs 
of overdensities, representing two richness clusters as those 
described above. The first is at redshift z=1.0 and the second 
at various higher redshifts. We found that, depending on the 
chosen threshold, the two structures can be separated if the dis- 
tance is greater then Az=0.15. However, once N galaxies are 
assigned to a single structure, the redshift of the latter can be 
determined with an accuracy of about Az/N 1 ^ 2 , which can be an 
order of magnitude smaller than Az. The above results refer to 
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the possibility of resolving (aligned) structures. It is worth not- 
ing that in real cases,once a (2+l)D density map is available, 
it is possible to adopt a multi-threshold technique, to identify 
possible physical substructures and/or projection effects. 

6. DISCUSSION 

It is clear that many methods exist for detecting clusters, 
depending on the type of the available data: from single-band 
images, where only brightness can be used to complement 
angular positions in the identification of overdensities like 
in matched filter methods (Postman et al. 119961 1. to images 
in a few bands where colours allow to identify the cluster 
red sequence (Gladders et al. 1998), to images in several 
bands allowing the determination of photometric redshifts, to 
spectroscopic redshift. Each method has his advantages and 
limitations. One of the main issues is our ability to compare 
the results of numerical simulations with galaxy catalogues 
obtained from different data sets. A general discussion of these 
problems has recently presented by Gal 2006 To establish 
whether one method is more or less efficient respect to the 
another in detecting structures is not a straightforward task. 
Given a mock catalogue, each method should be optimised, 
by a proper choice of the relevant parameters, to obtain a 
meaningful comparison. Moreover, mock catalogues obtained 
by different assumptions on cosmic evolution of the spectral 
and clustering properties of galaxies would lead to different 
optimisations. The problem of using the photometric redshifts 
to identify structures in the galaxy distribution has been 
tackled recently by Botzler et al. 2004 who discuss why stan- 
dard FOF methods, suited for spectroscopic redshift surveys, 
cannot be simply applied to photometric redshift surveys 
without a proper account of the large inherent uncertainties. 
They propose an extended EXT-FOF method which applies 
a two-dimensional FOF algorithm to slices of the galaxy 
catalogue, which are defined by photometric redshifts,taking 
into account their uncertainty. Clearly photometric redshifts 
add crucial information for the identification of structures 
in the galaxy distribution through the fitting of the observed 
spectral energy distribution to either theoretical or observed 
galaxy templates, evolving in cosmic time. Here we propose 
and discuss a different way of using photometric information: 
instead of comparing the distance between galaxy pairs, our 
method uses the statistical information about how many galaxy 
are in the neighbourhood of a given point and estimates a 
local density without introducing a fixed smoothing scale. A 
comparison of our method with EXT-FOF or other methods 
is beyond the aim of the present work and will be the subject 
of future investigations. What we want to stress here are only 
the new informations we obtain with our method: i) the spatial 
resolution is maximised in each point, since near density 
peaks the volume taken into account to evaluate the density 
is small, while it becomes large where the density is low; 
ii) once the density map is available, it can be analysed in 
different ways without re-running the algorithm. For instance 
the map can be sliced at different density thresholds to see 
how structures which appear separated at high density merge 
at lower densities. In principle it is also possible to apply an 



harmonic analysis to density maps. 

We have also obtained some new results by applying our 
(2+l)D method to real data in the Chandra Deep Field South 
which is one of the most studied fields in the sky. The existence 
of deep X-ray data from the Chandra observatory provided a 
sample of Active Galactic Nuclei (AGN) in the field (Giacconi 
et al. 2002). The spectroscopic follow up has shown the exis- 
tence of large-scale structure (Gilli et al. 2003). AGN, thanks 
to their intrinsic luminosity and to strong emission lines, are 
the ideal tracers of the mass distribution at the maximum dis- 
tance reachable by optical spectroscopy, under the assumption 
that their spatial distribution mimics that of normal galaxies. 
The same depth is either unreachable or requires huge expo- 
sure times for normal galaxies. Photometric redshifts, though 
with the limitation imposed by their poor accuracy, permit to 
identify overdensities of normal galaxies at the depth of AGN 
samples. This allows to study the relation between the space 
distribution of AGNs and galaxies. In the case of the CDFS, 
our photo-z based density detects two clusters, at z~0.7 and 
z~l, already identified in the spectroscopic redshift survey of 
Gilli et al. 120031 A third peak at z~ 1 .5 in the galaxy density 
distribution likely corresponds to a peak which is found in the 
spectroscopic redshift distribution of AGNs but not of galaxies. 
The analysis of deeper photometric data, which is in progress 
(Trevese et al.2006 ), is necessary to confirm the reality of this 
structure. Once galaxies belonging to a cluster are identified, it 
is possible to study the galaxy type distribution as a function of 
the local density. This provides an evidence of environmental 
effects on galaxy evolution and a way to quantify these effects. 
In the case of our analysis this allowed to prove, in a cluster of 
redshift 0.7, that the fraction of red galaxies increases with den- 
sity, as already known for lower redshift structures (Dressier 
1980 and Dressier et al. U997b . A similar effect appears at red- 
shift 1, though the evidence in this case is marginal. Clearly the 
density computed by our method is not the "real local density", 
due to the strong smoothing in the z direction caused by the 
redshift uncertainty and to the assumptions adopted to com- 
pensate the increasing loss of faint galaxies at high redshift. 
Thus any cosmological application of our method, devoted to 
study the evolution in cosmic time of galaxy clusters properties 
and environmental effects, will require a calibration through a 
detailed comparison with simulated catalogues. Obviously the 
same is true for any other cluster finding technique. However, 
there are properties of the cluster galaxy population which, to a 
first approximation, do not need a comparison with mock cat- 
alogues to provide significant physical information. This is the 
case of the slope of the cluster red sequence which has been 
traced up to z=1.27 (Blakeslee et al 2003 and refs. therein). 
According to Kodama et al. 1998 the evolution of the colour- 
magnitude (C-M) slope depends on relative age variations in 
early-type galaxies with different luminosities. Accordingly, 
the constancy of the slope up to redshift 1.27 is consistent with 
passive evolution of an old stellar population that was formed 
at high redshift. In this scenario, changes of the C-M relation 
slope are expected at redshifts approaching the star formation 
phase. An alternative interpretation of the constancy of the C- 
M relation, in the framework of hierarchical models of galaxy 
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formation, has been proposed by Kauffmann & Chariot 1998 
In this scenario metals, once formed, are more easily ejected 
from smaller disks. Large (bright) ellipticals are more metal- 
rich because they are formed from the mergers of large disks. 
In selecting rich clusters at high redshift one is biasing sam- 
ples towards objects that merged at the highest redshifts and 
for this reason they appear to follow the passive evolution. In 
both scenarios the C-M relation is due to a mass-metallicity re- 
lation and the evolution of the C-M slopes at high redshift con- 
tains critical informations on the origin of the C-M slope itself. 
Our analysis supports the evidence of a constant C-M slope up 
to z=l, meaning a correspondingly constant mass-metallicity 
relation. We stress that most of the highest redshift clusters 
detected so far were selected in the X-ray band, in particular 
those at z=1.24 (Blakeslee et al. 120031 Rosati et al. l2004l i and 
z=1.4 (Mullis et al. 120051 . However, both this clusters have a 
velocity dispersion of about 800 km s _1 and intracluster gas 
temperature kT ~ 6 Kev.typical of rich clusters (Bahcall 1988 
Arnaud et al. 2005). The X-ray luminosity of these two clus- 
ters in the 0.5-2 Kev band are 1 .9 x 10 44 and 3.0 x 10 44 erg s" 1 
h~l respectively (Rosati et al. 120041 Mullis et al. 120051 1. again 
typical of clusters with richness class greater than 2 (Ledlow et 
al. [2003 ). Rich clusters of these redshifts are just at the limit 
of our present (2+l)D analysis based on photometric observa- 
tions. Work is in progress (Trevese et al. 2006) to extend the 
study to a deeper samples (Kab ~ 27) which will allow the 
analysis at z~1.5 of even poor structures of the type of those 
detected by photometric redshifts in the present work in the 
CDFS. These less pronounced structures are hardly detectable 
in X-rays, due to the strong dependence on richness of the X- 
ray luminosity (Ledlow et al. 2003), and their C-M slope could 
differ from that of the richest clusters, if it constancy is only 
apparent and mainly produced by the cluster selection bias to- 
wards richer structures. 

7. CONCLUSIONS 

We have presented a new method to detect local overdensities 
in the galaxy distribution, based on: i) photometric redshifts, 
with proper account of the relevant uncertainty, to evaluate dis- 
tance; ii) distance to the n-th neighbour to evaluate densities. 
From a methodological point of view, we can conclude that: 

- although our (2+l)D method is limited by the redshift un- 
certainty, so that in principle only structures separated by 
more than 6 p i m , can be resolved, it still allows to dramati- 
cally increase our possibility of detecting structures respect 
to 2D analyses, as shown by the example in Figure[9] 

- if we restrict to structures of the type of Abell clusters, 
whose number density is about 10~ 5 Mpc~ 3 (Bahcall 1988), 
the average inter-cluster distance results d c \ ~100 Mpc , 
namely 8 p h i < d c i, meaning that the use of photometric 
redshifts allows to resolve even aligned structures, as shown 
by the simulation in section 5. 

From an astrophysical point of view, our analysis of deep 
multiband photometry in the CDFS field allows the following 
conclusions: 



- we have detected two structures at redshifts 0.7 and 1.0, 
whose existence was known from previous spectroscopic 
studies; 

- a third structure at redshift 1 .5 has also been detected but 
requires deeper data for confirmation; 

- the fraction of red galaxies inside the structure at z= 1 indi- 
cates a marginal density dependence while in the structure 
at z=0.7 the increase of the red fraction with density is seen 
very clearly; this extends the results of Dressier et al. 1 19971 
Carlberg et al.[2002and Tanaka et al. 120051 

- our results add new evidence in favour of constant slope 
of the C-M relation in clusters, at least up to z= 1 implying 
a constant mass-metallicity relation, according to the stan- 
dard interpretations; 

- our analysis shows that the average C-M relation for galax- 
ies belonging to the red population is consistent with a lin- 
ear extrapolation of the relation found by Bell et al. 2004 
up to z=1.5, complementing the results of Giallongo et 
al. 120051 who have shown that the blueing of the average 
colour of the red galaxies extends to z» 2 - 3; 

- the use of photometric redshifts will allow to analyse the 
redshift dependence of the C-M relation at high redshift, 
even in moderate overdensities, providing constraints on 
the very origin of the C-M relation. 

Thus, in spite of the low resolution in distance, (2+ 1)D analysis 
based on photometric redshift is an extremely valuable tool to 
complement other cluster finding techniques and perform large 
scale structure studies based on photometric surveys which, 
at present, posses unique capabilities in combining depth and 
field width. 
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